THE EFFECT OF PENGUINS IN THE Ba J/tpK^ CP ASYMMETRY 
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Performing a fit to tlie available experimental data, we quantify the effect of long-distance contri- 
butions from penguin contractions in B° — > J/ipK° decays. We estimate the deviation of the mea- 
sured iScp term of the time-dependent CP asymmetry from sin 2/3 induced by these contributions and 
by the penguin operators. We find A5 = <Scp(J/V'-K')-sin2/3 = 0.000±0.017 ([-0.035,0.033] @95% 
probability) , an uncertainty much larger than previous estimates and comparable to the present sys- 
tematic error quoted by the experiments at the B-factories. 



The measurement of the phase of the B^, — Bd 
mixing amplitude, given by twice the angle (3 
of the Unitarity Triangle (UT) in the Standard 
Model (SM), is one of the main successes of B- 
factories, and a crucial ingredient to test the SM 
and to look for new physics. The golden mode 
for this measurement is given by B'^ —> J/ipK'^ 
decays 0. These modes give a value of sin 2/3 
which is considered practically free of theoretical 
uncertainties and thus serves as a benchmark for 
indirect searches for new physics. Indeed, new 
physics can reveal itself by comparing different ob- 
servables - which all determine sin 2/3 in the SM - 
to the reference value from the J/tpK'^ modes. For 
instance, sin 2/3 can be extracted from the UT fit 
or from 6 — s- s penguin-dominated modes such as 
B — > (pKs or B ^ fl'Kg. Actually, possible hints 
of a discrepancy are being seen in both cases [1,1^. 

Impressive progress has been recently achieved 
at the S-factories in the measurement of the coef- 
ficient Sep of the time-dependent CP asymmetry 
in B^ J /ij^K^ decays. The experimental error 
on Sep has been pushed down to ±0.028 (statis- 
tical) ±0.020 (systematic) 0. On the theoretical 
side, previous estimates of the uncertainty in the 
extraction of sin 2/3 from iScp gave results below 
10"'^ (for a recent study, see ref. 0|) and therefore 
completely negligible. In this paper, we reana- 
lyze this issue with a new approach, described in 
detail below, obtaining a substantially larger un- 
certainty comparable to the present experimental 
systematic error. 

The decays of neutral B mesons into J/ il^K^ fi- 
nal states are dominated by a tree-level amplitude 
proportional to VcbV*^. Assuming the absence 
of additional contributions with different weak 
phases, it is possible to extract the value of sin 2/3 
from the coefficient Sep of the time-dependent 



CP asymmetry in these decays. As already men- 
tioned, the identification of iScp( J/V'^s/l) with 
sin 2/3 is affected by a theoretical uncertainty, com- 
ing from the presence of additional contributions 
having a different weak phase and possibly a rel- 
ative strong phase with respect to the dominant 
contribution j^]. Using the OPE, we write the 
expression of the decay amplitudes arranging all 
the contractions of effective operators into renor- 
malization group invariant parameters . In this 
way, we have 

A(5" ^ j/4,K'')=v:fyUE2 - P2) 

+v:,vup^'^ - P2) , (1) 

where E2 represents the dominant tree contribu- 
tion and the other terms are penguin corrections. 
Although three parameters (£'2, P2 and P^^^) 
enter the amplitude, for the purpose of this pa- 
per they can be treated as two effective parame- 
ters E2 - P2 and ^2^^"^ - P2. Neglecting the dou- 
bly Cabibbo-suppressed combination P^^^ — P2, 
a penguin pollution could come from P2. Even 
though this contribution might have an impact 
on the branching ratio, it certainly does not affect 
the CP asymmetry, since the two amplitudes carry 
the same weak phase. Conversely, because of the 
weak phase of Vub, P^™ — P2 might produce an 
effect on Sep and Ccp, although the impact on 
the branching ratio is expected to be very small. 

Being doubly Cabibbo suppressed, the value of 
^2^^^^ — P2 is hardly determined from B J/ipK 
decays alone. Therefore one needs to extract the 
range of this parameter from a different decay 
in order to study the impact of such a subdom- 
inant effect on sin 2/3. Indeed, the induced un- 
certainty on iScp increases with the upper bound 
of this range. It is then of the utmost impor- 
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tance to quantify this upper bound in a reliable 
way. Previous detailed discussions of the uncer- 
tainty AS = ScpiJ/ipK) — sin 2/3 have estimated 
the effect of P^^^ - P2 using the BSS mecha- 
nism recently supported by QCD factoriza- 
tion, to express penguin contractions in terms of 
local four-fermion operators 1^. However, QCD 
factorization holds only formally for this channel 
. Clearly, the importance of this measurement 
for testing the SM and looking for new physics 
calls for a more general assessment of the theoret- 
ical uncertainty. In the present work, we aim at 
providing a model-independent estimate of AS. 

To fulfill our task, we proceed in three steps: 
i) Neglecting Pj^™ ~ P2, we extract the absolute 
value of E2 ~ P2, using the experimental value 
of the branching ratio, ii) We extract [£'2 — ^2!, 
-^2^^^^ — P2 1 and the relative strong phase 5p from 
a fit to the SU(3)-related (up to the assumption 
discussed below) channel J/ tpir'^ . In this de- 

cay mode, P^^^ — P2 is not doubly Cabibbo sup- 
pressed and can be determined with good accu- 
racy. At the same time, we can compare the value 
of E2 — P2 obtained in the two channels to test the 
SU(3) invariance and the additional assumption. 
We can then take the range of Pg^^'^ — P2 from 
this fit (at 99.9% probability) as a reliable esti- 
mate of the range to be used in P" J/i/jK'^. iii) 
We repeat the first step, varying P^™ — P2 in the 
range obtained in the second step. In this way, we 
get the distribution of Sqp , to be compared with 
the input sin 2(3 to obtain AS. 

Let us provide some details about the second 
step. Using the same formalism of Eq. we can 
write the decay amplitude of P° J /ipir^ as: 

A(P° ^ J/7/'7r")=t/;t4<i(P2 - P2) 

+V:MP?'^~P2), (2) 

where all the combinations of CKM elements 
now are of the same order of magnitude and the 
additional (OZI-suppressed) contribution of the 
emission-annihilation EA2 parameter has been ig- 
nored .15.]. Even though the SU(3) symmetry is 
not exact (so that assuming the parameters to be 
the same in the two fits would require a difficult 
estimate of the associated error), we think that 
SU(3) is good enough to give us a reasonable es- 
timate of the allowed range of |P2^^'^ — ^2|- 

In the three fits, we use as input the determi- 
nation of the CKM matrix obtained by the UTfit 
Collaboration discarding the bound on p and f/ 



TABLE L Input values used in the analysis. All di- 
mensionful quantities are given in GeV. 





0.27 ±0.08 


pB^K ipB^^ 


1.2 ±0.1 


f.J/ii 


0.131 


ms 


5.2794 


P 


0.207 ± 0.038 


V 


0.341 ± 0.023 


A 


0.86 ± 0.04 


A 


0.2258 ±0.0014 


Gf 


1.166 ■ 10"^ 


(y,em 


1/129 



from po J/^K° 1^. To give a reference nor- 
malization factor for all the results, we use the 
value of P2, computed using naive factorization. 
All the inputs used in the fit are summarized in 
Tab. n We assume flat distributions for F^^^ 
and for pB~^K j pB^ir j^-^ ^y^q ranges specified [lfil |. 

Using the experimental value of BR(P° — > 
JJjhK^), we bound the absolute value of P2 — P2 
|1(tI |. Using the statistical method of UTfit 11], 
we assign a flat a-priori distribution to the abso- 
lute value |p2 — -F2 1 in a range large enough to fully 
include the region where the a-posteriori distribu- 
tion is non-vanishing. In this way, we reproduce 
the experimental value of the branching ratio with 
an indication of a significant effect of nonfactoriz- 
able corrections in IP2 — ^2!, as already noted in 
[13. We obtain IP2 - ^2! = 1.44 ± 0.05. Notice 
that, in the single-amplitude approximation used 
in this first step, the predicted Ccp is exactly van- 
ishing while iScp is, as expected, equal to the input 
value for sin 2/3 {Sqp = 0.729 ± 0.042). 

We now extract P^™^ - P2 from P° ^ J/ipir^. 
For this fit, we use the same approach but we 
retain in the amplitude jPa - -P2I, 1^2^™ - ^2! 
and the relative strong phase 5p. Together with 
the experimental information from the branching 
ratio and Cqp, we impose the constraint coming 
from5cp We allow |P2-^2| and jPp^'^-Pjl 
to vary in a range larger than the support of the 
output distributions, and 5p G tTjTt]. The re- 
sults are given in Tab. ^ As can be seen from 
the correlation plot in Fig. ^ two solutions are 
possible, with IP2 -^^2! and jPj^^'^ -P2I exchang- 
ing roles. Comparing the results of this fit with 
the value for |p2 — ^2! obtained from P — > J /-ipK^ , 
it is evident that only the first solution in Tab. Hll 
(with IP2 - -P2I = 1-22 ± 0.15) is compatible with 
SU (3) and with our expectations on the relative 
sizes of P2, P2 and P2™. Assuming therefore 
that this ambiguity is resolved in favour of the 
first solution, we repeated the fit with the cut 
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TABLE 11: Results of the fit of B° ^ J/i/)7r° (see the 
text for details). 



^CP 


0.09 ±0.19 


ricxp 
'^CP 


0.12 ± 0.24 


cth 
'-'CP 


-0.47 ±0.30 


'^CP 


-0.40 ± 0.33 


\E2-P2\ < 


1.22 ±0.15 




0.43 ± 0.43 


0.15 ±0.15 




2.87 ±0.43 


5p 


(-24 ±41)° 
(-146 ±50)° 








FIG. 1: Correlation between the hadronic parameters 
\E2-P2\ and \P^^^ - P2I, as obtained from the fit to 



jpGiM _ ^ 2\E2- P2\- The results are pre- 
sented in Fig. Inland in Tab. lIIII We underline the 
good agreement between this result and the deter- 
mination of |i?2 — ^2! from B — > J/^fjK'^, and we 
conclude that there is no evidence of large SU{3)- 
breaking effects or emission-annihilation contribu- 
tion in the determination of this parameter. We 
thus decide to use as input for the determination 
of AiS in B ^ J/^Ks a uniform distribution in 
the range [0, 2.3] for \P^™-P2\. This corresponds 
to the 99.9% probability range for \P^^^-P2 \ ob- 
tained in the fit. 

Repeating the fit of J/tpK" with the ad- 

ditional contribution of P^^^ — P2 in the range 
obtained above, we get the results in Tab. IIVI 
We also show in Fig. 13 the output p.d.f. for 
jpGiM _ g^jjj together with the difference 



TABLE III: Results of the fit of B" J/^tt° with the 
cut \P^™ _ P2 1 < 2|_E2 - P2 1 (see the text for details) . 

Ccp 0.09 ±0.19 CcpP 0.12 ±0.24 

Sep -0.58 ±0.24 Sep -0.40 ±0.33 

\E2-P2\ 1.22 ±0.15 |P2'^'^--P2| 0.38 ±0.38 

Sp (-34 ± 41)° U (-144 ± 19)° 



TABLE IV: Results of the fit of B° J/xjjK° (see 
the text for details), ^gp* (i^cp) represent the input 
(output) values of <Scp respectively. 



•-CP 


0.00 ± 0.02 


•^CP 


-0.01 ±0.04 


COUt 

'-'CP 


0.73 ± 0.05 


cin 
'-'CP 


0.73 ± 0.04 


1^2 -P2I 


1.44 ± 0.05 


|pGIM _ 


- ^^2 1 , 5p: see text 



A5. The result is 

A5 = 0.000±0.017 ([-0.035, 0.033] @ 95% prob.) . 

(3) 

Notice that, as anticipated, jPj^^'^ — ^2! and 
5p are poorly determined in this fit. In particu- 
lar. Fig. shows how the bound on the range of 
|pGiM _ p^i po j/^^^o extremely ef- 

fective in cutting out a long tail at large values of 
jpGiM _p^|^ thus reducing the uncertainty on AS*. 
Without this additional information, \P2™^ — P2I 
could have reached much larger values and corre- 
spondingly we would have obtained values of AS 
of order one. 

Had we boldly borrowed from the previous 
step not only the range but also the shape of 
jpGiM _ P2I, we would have constrained the de- 
viation of Sep from sin 2(3 even more, obtaining 
a value AS = 0.020 ± 0.007. However, given 
the theoretical uncertainties related to the SU (3) 
breaking and the neglected emission-annihilation 
contribution, this result is quoted for illustration 
only, and should not be used for phenomenology. 
A more reliable result can be obtained by adding 
a 100% error to the SU{3) relation between the 
hadronic parameters in the two channels. In this 
way we obtain 

AS = 0.001±0.015 ([-0.025, 0.026] @ 95% prob.) , 

fully compatible with our main result in eq. ((JJ. 
We conclude that our approach of extracting from 
B J/ip7T° the range of |P2^^^--P2| to be used in 
B J/tpK^ is fully consistent and does not siz- 
ably overestimate the error in AS. We also stress 
the importance of improving experimental results 
on B —> J/'07r'^ in order to reduce the uncertainty 
in the extraction of sin 2/3 from B J/tpK'^ de- 
cays. 

Our estimate of the error in AS is more than an 
order of magnitude larger than previous estimates 
and comparable to the present experimental sys- 
tematic error. This uncertainty should therefore 
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FIG. 2: Output distributions of hadronic parameters \E2 — P2\ (top left), |P2'^^~^2| (top right) and 5p (bottom), 
as obtained from the fit to B° J/ipiv^ with the cut \P^'^^ - P2I < 2\E2 - P2\ (see the text for details). 




IP" -P I 





0.05 0.1 

AS(J/vK°) 



FIG. 3: Output distributions of hadronic parameters P^™ (left), Spcim (middle) and A<S (right). 



be included in the value and error of sin 2/? ex- 
tracted from S^p. We believe that additional ex- 
perimental information on the decay modes con- 
sidered in our analysis will allow to reduce the 
uncertainty in AiS using the new method sketched 
in this paper and without any need of additional 
theoretical input. 
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